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This paper reviews the major sources of harmonic currents excited by an electromagnetic wave 
in an iostropic plasma. The discussion is based on the assumption that a plane wave of prescribed 
form interacts with a plasma whose response is described by the usual macroscopic equations for the 
electron velocity, density and temperature. Nonlinear interactions arising from (i) the magnetic field 
of the plane wave, (ii) the dependence of the collision frequency and other rate parameters on the 
electron temperature, and (iii) the oscillating space-charge induced by the wave in regions of inhomo- 
geneity, are examined in sufficient detail to provide an estimate of their relative importance in pro- 
ducing harmonic currents in various situations. The physical origins of these currents are discussed 
and a characteristic perturbation parameter provided for each of the interactions. 



1. Introduction 

When a high-frequency electromagnetic wave inter- 
acts with a plasma, there are various weak non- 
linearities through which harmonic currents can be 
excited, and these currents may, in turn, give rise 
to secondary radiation fields. Although certain re- 
lated nonlinear interactions have been studied in 
detail, mostly in connection with wave interaction 
effects in the ionosphere [Ginsburg and Gurevich, 
1960], the literature on harmonic excitation itself is 
relatively sparse and widely scattered. The reason 
for this is fairly obvious: the nonlinear sources re- 
sponsible for the appearance of harmonic currents 
are weakly excited and are therefore of little interest 
under most circumstances. Put another way, the 
linearized theory of wave propagation in a plasma 
works extremely well. However, just as the avail- 
ability of intense laser fields at optical frequencies 
has led to an interest in the nonlinear optics of crystal- 
line media [Armstrong et al., 1962; Bloembergen, 
1963], so also should the increasing power of radio 
and radar transmitters and the long path lengths in- 
volved in space communications lead to increased 
interest in the relatively weak excitation of harmonic 
signals in ionized media. 

Up to the present time the three main practical 
motives for studying the excitation of such signals 
have been: (1) the production of spurious signals in 
the ionosphere [Feinstein, 1950; Forsterling and 
Wuster, 1951; Vilenskii, 1953; Wetzel, 1963; Rydbeck, 
1963] ; (2) the design, or explanation, of useful harmonic 
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generators [Maiman, 1957; Uenohara et al., 1957; 
Baird and Coleman, 1961; Swan, 1961]: (3) the pos- 
sibility of diagnostic applications [Wetzel, 1961a; 
Rydbeck, 1963]. Harmonic spectra have been ob- 
served experimentally in a variety of circumstances, 
but their origins are seldom clear. For example, in 
those cases in which a steady magnetic field was 
present, exciting fields at the gyromagnetic frequency 
of the electrons have produced observable harmonic 
signals, but while the cavity measurements made by 
Maiman [1957] are consistent with a simple and 
reasonable theory, the more complicated computer 
calculations of Whitmer and Barrett [1961] are at 
variance with the experimental results in several re- 
spects. Harmonics produced in a microwave dis- 
charge without a steady magnetic field have been 
ascribed in one case to strongly inhomogeneous electric 
fields [Swan, 1961] and in another to modulation of 
the ionization frequency [Baird and Coleman, 1961]: 
in the latter case the proposed hypothesis accounts 
very well for the observed data, but in the form used 
it appears to lack any convincing theoretical basis. 
The general picture appears to be that while detect- 
able harmonics can be produced quite easily, the non- 
linear processes underlying their production are not 
usually well understood and most probably occur in 
combination. 

In an effort to bring some order to this picture, we 
will attempt in this paper to examine a certain class 
of nonlinearities which might be responsible for some 
of these harmonic spectra. Because of the com- 
plexity of the analysis in the anisotropic case, we will 
restrict our attention to isotropic plasmas. A steady 
magnetic field can be expected to enhance certain of 
the isotropic interactions in the neighborhood of 
cyclotron resonance, as well as to introduce additional 
harmonic terms having a different character. When 
appropriate, we will mention such effects in passing. 
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2. Basic Formulation 

Several authors have examined certain features of 
the various nonlinear interactions between a plasma 
and an electromagnetic field on the basis of kinetic 
theory, using the Boltzmann equation [Margenau 
and Hartman, 1948; Vilenskii, 1953; Fain, 1955; Rosen, 
1961]. This approach has definite virtue when study- 
ing a single interaction in detail, but is far too cumber- 
some if one wishes to evaluate the comparative effects 
of several nonlinearities within a given mathematical 
framework. Accordingly, we will choose a simpler, 
and physically more transparent, description of the 
plasma in terms of the particle and momentum trans- 
port equations, augmented by a phenomenological 
energy equation. The primary electromagnetic wave 
will be assumed, when possible, to have the form of 
a transverse plane wave with prescribed local 
parameters. 

The total current density due to charge motion in 
the plasma is defined as 



j = ~e(n e Ve — riiVi) 



(2.1) 



where — e is the electron charge and n e , v e and rn, v\ 
denote the density and average velocity of the electrons 
and ions respectively. Since the densities are roughly 
the same and the ion velocity is smaller than the elec- 
tron velocity by the factor m/M, the mass ratio of elec- 
tron to ion, we may neglect the ion current and regard 
the ions as merely providing a uniform neutralizing 
background. Henceforth, then, we will be interested 
only in the electron current which, dropping sub- 
scripts, will be written 



= — env. 



(2.2) 



Harmonic components in this current will arise from 
harmonic components in n and v , so we must examine 
the nonlinearities in the equations governing the 
response of the plasma to the applied electromagnetic 
field; namely, the equation for the electron velocity 



- — \-v Vi; = (E + vXB) — Vc(T e )v; 

dt m mn 



the continuity equation 



dt 



(2.3) 



(2.4) 



and, to close the system, an equation for the electron 
energy (written here as a phenomenological equation 
for the electron temperature [see, e.g., Ginsburg 
and Gurevich, I960]): 



dT e 
dt 



2e_ 
3k 



vE-&p c (Te)[Te-To\. (2.5) 



It has been assumed in the moment equation (2.3) 
that the electron velocity distribution is isotropic, 



so the divergence of the stress tensor has been written 
as the gradient of the isotropic pressure p ( > = nkT e , 
where k is the Boltzmann constant. In the last term 
of (2.3) the temperature dependence of the effective 
collision frequency is noted explicitly. The tempera- 
ture equation (2.5), if multiplied by (3/2)nk, simply 
states that the rate of increase of the average electron 
energy per unit volume is equal to the difference 
between the electric power delivered to the electrons 
and the rate at which energy is lost by collisions with 
heavy particles held at the equilibrium temperature 
To, 8 being the energy loss parameter (d= z 2m/M for 
elastic collisions, where M is the mass of the heavy 
particle). 

Rather than expand everything in Fourier series 
and then attempt to disentangle the sources of har- 
monics from the resulting confusion, we will seek 
instead to identify the individual sources of harmonics 
in physical terms and to provide a characteristic 
measure of the relative importance of each. These 
sources can be anticipated from an examination of 
the equations (2.3) to (2.5). The plasma is driven by 
the prescribed field E (and B) and its response appears 
in the variables v, n, and T ( >. We know even from 
simple orbit theory that the magnetic field B of the 
plane wave could perturb the trajectory of the elec- 
tron in such a way that its motion relative to the 
phase plane of the wave might be anharmonic. Such 
a nonlinearity would be "intrinsic" in the sense that 
it would always be present when a wave propagates 
in a plasma, and would thus provide a threshold 
of interest in evaluating the importance of other 
nonlinearities. 

A nonlinearity like the v ■ V v term in (2.3) cannot 
itself be regarded as a primary source of harmonics, 
although it will affect the detailed response of the 
plasma to those primary sources which produce 
changes of electron velocity along the streamlines. 
Behavior of this latter type might arise from pressure 
gradients produced by a harmonic temperature wave 
traveling with the primary wave, or from space-charge 
variations induced when the primary wave encoun- 
ters an inhomogeneity. Moreover, the existence of 
temperature oscillations could introduce harmonic 
components into v through the 7Vdependence of the 
collision frequency in the last term of (2.3). 

Our general approach will be to use the formula- 
tion provided by (2.3) to (2.5) to find a semiquantitative 
description of the effect of these various sources, treat- 
ing each as if all the others could be ignored. It will 
be convenient to have a zeroth-order approximation at 
hand, representing the linearized solution to (2.3) 
to (2.5) for the case of constant electron temperature, 
pressure and density. If the primary transverse 
plane wave is written as 



~S(~r\ . t) =£*?-«"* cos (at — /3k • "?) 



(2.6) 



where £0 is a constant vector, k is the direction of 
propagation {E • £ = 0), oj is the wave frequency and 
fi and a are the propagation and attenuation constants 
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in a uniform plasma of density n {) , then the solution of 
(2.3) in the zeroth-order approximation becomes 

-*n» (e/m)Eoe- a "? 



(to 2 + i/* 



[w sin (cot — (3k • r) + v c cos (cot — /3k • r )] (2.7) 

and from (2.2) the corresponding induced current is 
j (0) = — eno v (0) . This zeroth-order current is just the 
linear polarization current produced by the primary 
wave in the plasma and thus determines the values of 
/3 and a for this wave. In addition to generating har- 
monic currents, the nonlinearities will also perturb the 
propagation of the primary wave by producing slight 
changes in /3 and a. 



and combine (3.1) and (3.2) into a single equation 



du 

dt 



l 6 

v v -\ E {) (z) cos (cot — Bz — 6) 

c p m 



u — E(z, t) 

m 



(3.3) 



for the complex variable a. (For ease of writing, we 
have used the abbreviation E (z) for E exp (— az).) 
By the use of an integrating factor, the steady-state 
(particular) solution of this equation becomes 



exp 



. eE (z) 
— v c t — i sin (cot — fiz — 4>) 



x^E { 

m 



*> / 



cos (cot — (3z) exp 



v c t — i 



.eE (z) 



3. Intrinsic Electric-Magnetic Field 
Interaction 

A free electron in the field of a plane wave experi- 
ences both electric and magnetic forces, and although 
these forces vary with time, it _is in the nature of a 
plane wave that the average . of ExB always points in 
the direction of propagation /c, so the electron should 
experience a drift in the K-direction. Since the trajec- 
tory of the electron is now curved in the plane of 
polarization (established by E and /c), the projections 
of its velocity on the phase plane will no longer be 
sinusoidal — in fact, on the basis of symmetry we would 
expect that a Fourier representation of the motion 
would contain odd harmonics only. 

These conclusions can be made more precise by 
solving a reduced version of (2.3) as completely as 
possible. Let a plane wave whose electric vector is 
given by (2.6) propagate in the z-direction with E {) 
= E {) e. r , so that E(z, t) = e.,E {) exp {—az) cos (cot — fiz). 
In rationalized MKS units, the associated magnetic 
field will have the form 

B(z, t) = e y c~ l E exp (—az) cos (cot — pz — <f>) 

where c p is the wave velocity in the plasma and the 
phase angle cf> takes account of the phase shift of the 
magnetic field in a lossy medium and is defined by 
tan (/> = a/p. 

An approximation to v which accounts only for the 

vxB effect can be found by neglecting v • V v and 
V p ( > in (2.3), assuming that v c is a constant, and writing 
the resulting reduced equation in coordinate form: 



dv. r 
dt 



+ v c Vz= £(z, f) + — v z B(z, t) (3 l) 

m m 



Under steady-state conditions there will be no v y com- 
ponent, so we may define a new variable u = v x + iv z 



sin (ojt — /3z — (f>) 



dt. (3.4) 



The exponential factors containing sin (ojt — /3z — (/)) 
can be replaced by the Bessel function expansion 



exp 



■ eEpiz) . " 

:/ sin ((x)t — Bz — (b) 

Cnmco 



= ^Jn(volc p ) exp \ ± iniojt - [3z - $)] (3.5) 

where J„ is the Bessel function of the first kind of order 
ai, and we define the parameter Vo = eEo(z)lmco. 

Upon performing the integration term by term, mul- 
tiplying the series together and separating the harmonic* 
terms into real and imaginary parts, the orthogonal 
velocity components v, r and v z may be written in the 
following approximate form: 

v. r = —V(){Jl\s\\\ (o)t — pz) + A cos (cot— fiz)\ 

+ (3^2-2^) sin 3 ( w/ -/^ -§</>) 

+ (-gJoJz-^JU A cos 3 L* -/&--<£ 



+ 



. (n odd) 



(3.6) 



v, = -\-v i) {j {) J l \(/\- i -\) sin $ + cos (/>] 
+ hJoJdcos 2((ot-pz—%tf))-A sin 2(tot- (3z- h(j>)] 



+ 



(n even) 



}. (3.7) 



Only the lowest order Bessel functions have been 
retained, and the solution contains only terms through 
first order in the parameter A = iv/co, which will hence- 
forth be assumed to be small. Note that the trans- 
verse velocity v x contains only odd harmonics, while 
the longitudinal velocity v z contains only even har- 
monics. It should be emphasized that these are 
solutions of (3.1) and (3.2), not of (2.3). The existence 
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of v z means that both v V v and V p e should be in- 
cluded in any attempt to solve (2.3). On the other 
hand, one expects that a perturbation solution of (2.3) 
based on the velocities given above would not alter 
the results in a qualitative way, leading at most to a 
correction of some of the numerical coefficients. 

Since the parameter Vo is the magnitude of the un- 
perturbed electron velocity in the case v c = Q, and c p 
is always greater than the speed of light in vacuum, 
the ratio Volc p is small for nonrelativistic velocities and 
we can use the small-argument approximations for the 
Bessel functions: 



by 



jxa = — e(noVxa + n 2 v x i) 



(3.13) 



with v x \ and v X 3 taken from the first two terms of (3.9) 
and n 2 taken from (3.12), all of them in the limit A = 0. 



JxZ 



24 



e 2 sin 3{cot — (3z) 



en () v ( ) 



sin {cot — /3z) cos 2{cot — fiz) 



Jn(Volc P ) = (j ^AVolc P )<<l. 



(3.8) 



12 



€ 2 (enoVo) sin 3{cot — (3z) 



If these approximate forms are used in (3.6) and (3.7), 
the expansions can be expressed as power series in 
the small parameter e = (volc p ): 



V.r ■ 



■■ — i>o[sin {cot — /3z) + A cos {cot — (3z)] 



~2i eVo 



sin 3{cot — /3z — — cf>) 



-h— A cos 3{wt-pz--cf)) 

*ei; [cos cf)-\-{A-* — A) sin </>] 



+ 0(e 4 ) (3-9) 



+ 7^o cos 2 loot — /3z — -cf) J 

-|Asin2(co*-j3z-!<M + 0(e : <). (3.10) 

(It should be remarked that sin c/>/A remains finite for 
A = 0.) However, according to (2.4) the appearance 
of a second-harmonic velocity in the longitudinal di- 
rection implies the excitation of a corresponding space- 
charge wave: 



dn> 



dv z2 



dt dz 



(3.11) 



where numerical subscripts are introduced to identify 
harmonic components of the various plasma 
parameters. 

Substituting from the second term of (3.10), we 
obtain 



+ - € 2 {enoVo) sin {cot 
o 



-Pz). (3.14) 



The result dramatizes the difficulties encountered in 
nonlinear analysis of this type, even when pursued from 
the naive point of view we have adopted here. The 
products of harmonic components contain both sum 
and difference terms, the latter feeding corrections 
back into the lower order components. Fortunately, 
however, if these corrections are viewed in terms of a 
characteristic perturbation parameter, they will 
always be of higher order than the dominant contri- 
bution to a given harmonic and should therefore be of 
little practical significance. 

The intrinsic electric-magnetic field interaction can 
therefore be expected to generate a harmonic spec- 
trum of currents whose dominant terms are given by 

jo =*- e{en v )[cos ^-HA -1 — A) sin cb] e z 
j\ — {erioVo) [sin (cot — /3z) -f A cos {cot — (3z)]e. r 



ji - 7 €(en v ) [cos 2 f cot — (3z — - cf) ) 



— -A sin 2 [cot 



-*-£»)} 



i _., 



= - — € 2 {en Vo) [sin 3(cot - f3z) + 0( A )] e. r 



J2n ^ e 2 ""^, 



n 2 



2 



cot-(3z--cf) 



— -A sin 2 [cot — fiz- 



¥) 



(3.12) 



The definition of harmonic currents will therefore in- 
volve not only the harmonic velocity components, but 
harmonic density components as well. To illustrate 
the implication of this, let us evaluate the third- 
harmonic component of the transverse current defined 



j2n + \ 



OC *:2n ( 



€*"ej 



(3.15) 



The third harmonic has been simplified in accordance 
with (3.14) and the last two entries summarize the 
strength and polarization of the even and odd har- 
monics respectively. The anticipated longitudinal 
drift current shows up in 70, while j\ is just the linear 
polarization current in the field of the primary wave. 
As expected, the transverse current contains only 
odd harmonics, and these decrease by a factor e 2 from 
one term to the next. These transverse currents can 
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generate odd harmonic waves traveling in the direc- 
tion of the primary wave, provided that the dispersion 
at the primary frequency is so small that /3.i = 3/3i 
over an appreciable distance. On the other hand, it is 
easy to show from Maxwell's equations that there is 
no magnetic field associated with the longitudinal 
currents, so there will be no radiated harmonic fields 
of even order. 

The situation is quite different in the presence of a 
steady magnetic field. If, for example, the direction 
of the steady field coincides with the magnetic field 
of the wave, a traveling space-charge wave of order e 
is produced at the primary frequency a> which, in turn, 
provides a source for a second-harmonic transverse 
wave. 

The procedure used to find v x and v z in this section 
resembles that used by Maiman [1957] for calculating 
harmonic currents produced by prescribed cavity 
fields. Moreover, some currents of this type were 
obtained by Chen [1962] for the case in which *v = 0. 
Chen's results are incorrect, however, containing a 
spurious factor (1 + cofjaj 2 ) which resulted, effectively, 
from the introduction of an internal polarization field 
in addition to the impressed field. 



4. Electron Temperature Effects 

There are many interesting and important phe- 
nomena associated with the dependence of the 
effective collision frequency of the electrons in a 
plasma upon the electron temperature T ( , but only 
a few authors have considered the role of this non- 
linearity in the production of harmonic currents 
[Vilenskii, 1953; Rosen, 1961; Visvanathan, 1962]. 
There are actually several other ways in which electro- 
magnetic field-induced electron temperature varia- 
tions can produce harmonic currents in a plasma, but 
since the modulation of the collision frequency is 
probably the most important, we will consider it in 
greatest detail. 

An approximate solution of (2.5) can be obtained by 
assuming that in this equation, at least, the collision 
frequency v c is constant. Under this condition, and 
with v replaced by its zeroth-order approximation 
(2.7), the solution has been shown by Ginzburg and 
Gurevich [1960] to take the form 



&(v c l<x>) << 1 and (4.1) simplifies to 

~? = 1 +^[1 - \8 cos 2(G>t-/3z)+ %8A sin 2(ajt-pz)]. 

(4.3) 

Now, the definition of the effective collision frequency 
v c involves the electron velocity distribution function 
(hence some average energy parameter which we take 
here to be the electron kinetic temperature T e ) and the 
force law governing the scattering interaction between 
the electrons and the heavy particles. Accordingly, 
it is convenient to adopt a heuristic expression for v c 
of the form 



vJTe)=V C0 (T e IT*Y 



(4.4) 



where v co is the equilibrium value v c {Tq). 

The parameter cr depends upon the nature of the 
scattering interaction. For example, cr= 1/2 for rigid- 
sphere collisions appropriate to low energy scattering 
by molecules, cr = for the induced polarization inter- 
action usually occurring at high energies, and cr=— 3/2 
for coulomb scattering by ions. In general, a value 
of cr for a particular molecule and T e can be estimated 
from empirical curves of collision probability versus 
electron energy. Many such curves are given by 
Brown [1959]. 

We now assume that Te/T >> 1, which allows us 
to ignore the first term on the right of (4.3) and write 
(4.4) in the form 

MT C )^^{1- i8[cos 2(a>*-£z)-A sin 2(ojt-f3z)]}' r 

(4.5) 



where 



v^=v cl ,(T E jT a r. 



(4.6) 



Since 8 << 1, equation (4.5) can be expressed as a 
power series in 8: 



,o" , 



v e (T e ) - 4 01 - 1/<» g 8 [cos 2(a>t - j3z) - A sin 2(u>t - fiz)] 

,o-(o--l) 



+ l/0> - 



T e = T +T E 



. , (8v* - 2w 2 )3 

1 + W+8^ COS2M "^ 



[cos 2{a)t-pz)-A sin 2(<ot - /3z)] 2 + . 



(4.7) 



yhere 



doiv(2 + 8)8 0/ , Q ,-, (A ,v 
--———-sin 2M-0z)] (4.1) 
(zco)- + o l v l c 



T E 



e 2 Ej{z) 



3/cmoW + v*) 



(4.2) 



If the temperature relaxation time r=l/8^ r is much 
larger than the period 27r/co of the exciting wave, then 



Mr e )-i>< 0) +8i/i>+8v;> + . . . 



(4.8) 



where the v c n) are identified by comparison with the 
coefficients of 8" in (4.7). 

Let us again take a reduced form of (2.3) in which 
v ' V v , V p e , and now v X B are all ignored, and attempt 
a perturbation analysis based on an expansion of v in 
the form 



v = v {0) + dv {l) + d 2 v (2) + . 



(4.9) 
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If we substitute (4.8) and (4.9) into 



— = E(z, t)- 

at m 



MTe)^ 



(4.10) 



and equate successive powers of 8, the coefficients in 
(4.9) are found to satisfy the equations 



3 (0) —> e _» 

— + ^ ( > (0) = -E(z, t) 

dt m 



dt 



— +v«>W 1) =-vi 1 ¥ 0) 



(4.11a) 



(4.11b) 



out, for example, that in the case of collisions with ions 
(o" = — 3/2) the currents can be obtained from (4.15)- 
(4.17) simply by replacing 8 and 8 2 with — 38 and— 158 2 
respectively. Incidentally, the perturbation of order 8 
shown in the primary current in (4.15) was found by 
Taylor [1961] and by Epstein [1960] using a more 
complicated formulation of the problem. 

Returning to (4.3), we see that a second-harmonic 
temperature wave accompanies the primary wave 
and ought to produce a longitudinal diffusion cur- 
rent through the Vp c > term in (2.3). The appropriate 
velocity components can be estimated by neglecting 
v ■ V v, v X B and v c (T e )v and substituting V p e — kn^7T e 
with T e given by (4.3). 

The solution is given by 



dv^ 
dt 



+ v!?¥ 2) 



,(1U1)_ „(2U0) 



■CT 



(4.11c) 



The solution for the unperturbed velocity v {0) is just 
(2.7) with K = e z and v c replaced by i/ 0) , while the solu- 
tion for each higher-order coefficient requires knowl- 
edge only of the lower-order coefficients. Sacrificing 
generality for clarity, we write below the solutions 
for the first three coefficients in the special case cr = •>: 

fo[sin (cot- 



i*o) . 



i^oA 



-fa) + ^ cos (cot -fa)] 

1 



(4.12) 



cos (cot — fa) — — cos 3(cot — fa) 



(4.13) 



^ + -t, A 



cos (cot — fa) — - cos 3(cot — fa 



+ — cos 5(cot — fiz) 



. (4.14) 



Terms of order A 2 are again ignored, and V{)=eE {) (z)lmco. 
The velocity coefficients are all transverse, so there 
are no oscillating density components and the total 
current is therefor simply j = — en v, which separates 
into the following harmonic components: 

jd -erwo[sin (cot -fa) + A(l + 8/8~8 2 /64) 

cos (cot- fa)} (4.15) 

; C 3 * - e/ioVo(8A/24) (1-8/16) cos 3(cot-fa) (4.16) 

j C 5 — enm(& 2 A/640) cos S(cot —fa) (4.17) 



jc(2n+l) ^ l^o8 W A. 



(4.18) 



Modulation of i^ c by the primary wave thus produces 
an odd harmonic spectrum of purely transverse cur- 
rents whose magnitudes fall off by a factor 8 from one 
to the next. An additional numerical factor enters 
from the series (4.7) and the subsequent integrations. 
This factor depends on the value of cr, and it turns 



Vz ^llK »otan<*,-i C t;o 



cos 2(cot — fa — -</)) 



3 1 

- ^ A sin 2(cot — fa — -cf>) 



(4.19) 



The second harmonic term is exactly of the same form 
as that found in the last section for the intrinsic inter- 
action. Higher-order longitudinal harmonics could be 
found by reinstating the v • V v term in (2.3) and pro- 
ceeding with a perturbation calculation. There is 
no point in doing it here. It is sufficient to observe 
that such terms would, as before, generate even har- 
monics of order e 2 in the density n as well as numerical 
corrections in the higher-order components of v 9 
but would not be expected to alter the general features 
of the transverse currents in any important way. The 
zeroth-order term in (4.19) is a diffusion velocity result- 
ing from the decay of the field amplitude (hence T e ) 
as the wave moves through the slightly lossy plasma. 
Although its existence is of interest, its form in (4.19) 
cannot be taken seriously. If the electrons are urged 
to flow in the direction of wave propagation by a temp- 
erature gradient, the resulting charge separation in- 
troduces space-charge forces and the problem becomes 
one of ambipolar diffusion. 

Finally, we should mention the effect of temperature 
oscillations on such chemical rate processes as ioni- 
zation, recombination, and attachment. These 
processes can be included by introducing energy 
dependent coefficients in rate terms added to the right 
side of (2.4): 



dn 
It 



+ V • (n v ) = Pi(T e )n - a r (T e )n 2 - v n (T e )n. 



(4.20) 



The rate coefficients are the ionization frequency 
vu the recombination coefficient a r and the attachment 
frequency v a - The perturbation procedure used 
earlier to find the effects due to v c {T e ) can be applied 
to (4.20) as well, and the results are much the same. 
As an example, we have calculated the dominant third- 
harmonic current due to the temperature dependence 
of the recombination coefficient a r . The result is 

jrz = C(en vo)d(v ro lco) cos 3(cot — fa) (4.21) 
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where C is a numerical constant and v r o is the zeroth- 
order recombination frequency defined by v,o 
= a r {T h )n {) . The other coefficients in (4.20) would 
contribute harmonic currents of similar form. 

The experimental results obtained by Baird and 
Coleman [1961] remain a puzzle. They found both 
odd and even harmonics radiated from a microwave 
discharge between two posts in a waveguide when an 
additional d-c electric bias field was applied. Their 
measurements were remarkably consistent with the 
predictions of a simple theory based on the assumption 
that Vi ~y\v\ on the right side of (4.20), where y is 
a constant and v is simply the zeroth-order electron 
velocity in the combined a-c and d-c electric fields. 
Unfortunately, there is no convincing justification for 
this assumption, (their simple argument using a one- 
dimensional velocity distribution function contains 
an error). Our formulation would also lead to both 
odd and even harmonics in combined a-c and d-c 
electric fields, but in order to obtain anything resem- 
bling the relation v\~ |v|, the condition on the solution 
of (2.5) would have to be hv v > o>, which is difficult to 
satisfy at microwave frequencies. 



5. Effect of Plasma Inhomogeneity 

Up until now, our discussion of nonlinear interactions 
has supposed that the plasma was homogeneous — 
at least to first order. If the plasma is inhomogeneous 
to start with, it becomes impossible to prescribe the 
primary wave without considering a specific case lor 
which a solution, generally of a complicated propaga- 
tion problem, is available. However, the essential 
feature of an inhomogeneous plasma is its nonvanishing 
density gradient, so we can get some idea of the char- 
acter of harmonic currents in plasma inhomogeneities 
by studying what happens when a quasi-plane wave 
encounters a region in which V/io ^ 0. 

In order to maintain a connection with our previous 
discussion, we suppose that the wave propagating in 
the inhomogeneous plasma may be written in the form 



E{r,t) = Eq(t) cos [cot — k(r) • r) 



(5.1) 



where k{r) is the wave vector. Such an expression for 
E is generally valid only in the limit of geometrical 
optics. The average velocity of the electrons under 
the influence of this field will be taken to be the zeroth- 
order approximation (2.7) in the limit of small A. 
To avoid confusion with v (0} in previous sections, we 
will write this as 



v — E {) {r) [sin {ojt — k • r) 

mix) 



+ A cos (cot-k- r)}. (5.2) 

Since V n () # we will need (2.4) at the outset, which 
requires that 



— +V/i-t' + /i W=0. 
dt 



In general, the second and third terms in (5.3) will 
not vanish and there must therefore be a time-varying 
component of electron density. Considering this com- 
ponent to be a perturbation on the original electron 
density n {) (r) of the inhomogeneous plasma, we write 



n = n {) (r) + n {l) (r,t) 
and substitute into (5.3), yielding 



dt 



- + Vrc 'l> , + /ioV-T' = 



(5.4) 



(5.5) 



as the equation for n {l) . Since it has been assumed 
that the original electron distribution rio(r) is perfectly 
neutralized by a corresponding ion distribution and, 
moreover, that the ions remain unmoved by the held, 
then n {l) must represent an oscillating space charge. 
It is easy to see how this space charge comes about. 
Let us count electrons in a small volume located within 
the inhomogeneity, keeping in mind that the ion popu- 
lation in this volume will remain stationary. If the 
primary field forces electrons to move down the density 
gradient, the volume will acquire a net negative charge 
as electrons arrive in it from regions of higher density. 
Conversely, if the electrons are forced up the gradient 
the electron population in the little volume comes from 
regions of lower density, and since the ions remain 
fixed, the net charge is positive. For a wave-like 
driving field there will be a wavelike space charge 
perturbation moving along the inhomogeneity. 

The existence of space charge implies the existence 
of a longitudinal component of the electric field, that 

is E must satisfy 



V-£(r7 0=-™ (1) /e () . 



(5.6) 



Since it is assumed that (5.1) is a proper solution of 
the propagation problem, we may use (5.6) and (5.2) 
to deduce that 



V-7=- 



€ {) mco z 



dn {l 



v v n K ' 



dt 



which, when substituted into (5.5), yields 



(5.7) 



(5.8) 



(5.3) 



where we have introduced the plasma frequency 
(ol = e 2 nole m. The steady-state solution of (5.8) when 
(Op < a) is found to be 

(i) (e/jn)V/io - Ep rr 

» (1> =V_ w p 2+KA)2 {[<*-«5>-«&*] 

cos ((ut — k-r)— oj 2 A sin (cot — k • r) } (5.9) 

or, if terms of order A 2 are ignored and the abbreviation 
e = ( 1 — w 2 /g> 2 ) is used, 



758-H50 0-65— 10 
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nU>s* 



eVrio ' E 



cos (a)t — k • r) 

A 



sin (cot — k - r) 



(5.10) 



It should be remembered that n , E , e and A; are all 
functions of position r. 

A nonlinear current is produced when these space 
charges interact with the field. Using (5.2) and (5.10) 
in j (1) = — en {1) v\ we get 



-£> 1 A wJ e 3 Vtt • £p 7f r* 1 e 3 Vno^Eo -> 



or m'co^e' 



(r) 



sin 2(o)t — k-r)-\ — A cos 2(<at — k • r) 



(5.11) 



The first term is a_d-c drift current in the direction 
of the field vector Eq\ the second term is^a second 
harmonic current, also in the direction of Eq. Notice 
that if the inhomogeneous region of the plasma con- 
tains a point of plasma resonance (to p = co), the current 
increases sharply in the neighborhood of this point 
due to the vanishing of e in the denominator. Under 
these conditions one must return to the more complete 
expression for n (l) given in (5.9). The second har- 
monic current in (5.11) is then large but bounded by 
collisions when co p = co, having the expression 



/(I) : 



1 e 3 Wo • E 

2 m 2 to 3 A 



Eq ( r) [cos 2 (cot — k-r) 
- 2 A sin 2 (tot - A*-~r) ] . (5. 12) 



In order to have any confidence in the evaluation of 
this type of nonlinear interaction in inhomogeneous 
plasmas for which to p ~ co, it is really necessary to 
solve (2.3) to (2.5), or some equivalent system, along 
with Maxwell's equations for the field. In the vicinity 
q£ plasma resonance, for example, the field amplitude 
Eo(r) will itself be strongly inhomogeneous, and there- 
fore the electrons will be subjected to a force of 
changing magnitude along their trajectories. This 
leads to strongly anharmonic oscillations over a 
narrow region and could thereby produce harmonic 
waves [Forsterling and Wiister, 1951]. The sinusoidal 
motion of an electron will be distorted only to the 
extent that its displacement from equilibrium is large 
enough for the electrons to "see" the field inhomo- 
geneity. 

If we indroduce a length A/, to measure the distance 
over which the amplitude of the electric field changes 
appreciably, and use (f) to denote the maximum 
displacement of an electron in the average field in 
the neighborhood of a point, then the ratio (£)/A# is 
a measure of the importance of this nonlinearity. 
Field inhomogeneity has occasionally been proposed 
as a source of harmonics in microwave discharges 
[Swan, 1961 and references cited there], but to con- 
sider it further would carry us beyond the scope of 
our present discussion. 



If the plasma is underdense, in the sense that co* 
is always much smaller than w 2 , then e — 1 and the 
amplitude of the primary wave varies slowly with 
position. This means that (5.2) can be integrated 
locally to give the electron displacement tj = fv'dt 
at 7? When A = 0, the magnitude of this displacement 
is just 



€ = 2* 



(5.13) 



Considering the case of strongest interaction in 

which Eq is parallel to V/io, and estimating the local 
nonuniformity of the plasma by a characteristic length 
A n defined by | Vn | — no/A n , the magnitude of the 
second harmonic current in (5.11) may be written as 



I i* 1 ) I 

\J2 I 



en 



eE 
mco 



eE Q 



2mco 2 A n 
enov')(£l2A n ) 



= |y 1 |(f/2A„) (5.14) 



where \ji \ is the magnitude of the primary current. 
Clearly, the perturbation parameter that emerges here 
is (£/2A n ), which is the ratio of the maximum elec- 
tron displacement in the field to the scale of the density 
gradient. This parameter is usually quite small. 
As pointed out elsewhere by one of us [Wetzel, 1963], 
the dominant harmonics produced by this interaction 
form a sequence in which the pth-harmonic current 
is j p = — en p -\v\ where the electron density component 
n p -\ is obtained from 



dn, 



P-l 



dt 



+ V-7 P ( i ) 1 = 0. 



(5.15) 



Under special conditions this sequence contains both 
odd and even harmonics which fall off from one term 
to the next by the factor (£/2A n ). It is clear, however, 
that the spectrum is determined by the way in which 
zio(r) changes with position along the direction of the 
primary electric field vector, so it is impossible to make 
any general statements about the relative amplitudes 
of the harmonic components. Various aspects of this 
type of nonlinear interaction in an inhomogeneous 
plasma, including the effect of a steady magnetic field 
on n {1 \ may be found in papers by Ginsburg [1958] 
and by Wetzel [1961a, b, 1963]. 

There is still another nonlinear current which has 
nothing to do with space charge, but is associated with 
inhomogeneities and properly belongs among the 
electron-temperature effects discussed in the last 
section. In an inhomogeneous plasma, the pressure 
gradient in (2.3) consists of two terms 



V (nJcTe) = n k V T e + kT e V n . 



(5.16) 



The first term was treated in section 4, but the second 
term gives rise to another second harmonic of electron 
velocity which, when multiplied by — en , provides a 
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current 

12 war 



sin 2(o>/ — k • r) 



+ « A cos 2(cL>t — k • r) 



(5.17) 



This current is the same size as the corresponding 
space charge current, but is directed along the density 
gradient rather than along the electric field vector. 
Its origin in (5.16) shows that it is a harmonic com- 
ponent of the diffusion current due to concentration 
gradients and should therefore be regarded with sus- 
picion, in view of the discussion following (4.19) 
concerning diffusion currents. 



6. Summary 

This review of the major sources of induced har- 
monic currents in an isotropic plasma has shown that 
these sources are quite weak and that the harmonic 
currents are all measured by small perturbation para- 
meters of sonic kind. We summarize the main results 
of the preceding sections below. 

A. Intrinsic E — B interaction: Odd harmonic trans- 
verse currents of relative magnitudes 1 7211+1 1 ^€ 2 "|ji|> 
where e = {volc p ) is the ratio of the unperturbed elec- 
tron velocity in the field of the wave to the phase 
velocity of the wave; even harmonic nonradiating 
longitudinal currents plus a drift current. The para- 
meter e can be estimated from the formula e — 10~ 4 
(Eolf), where Eq is given in volts/meter and / is given 
in megacycles. 

B. Electron temperature effects: Odd harmonic 
transverse currents of relative magnitudes |j2n+-i | 
0C 8 w A|7"i I , where 8 is the energy loss parameter and 
A is the ratio (*v/o>); even harmonic longitudinal cur- 
rents resembling those under A. above. It should be 
noted that the results of section 4 are based on the 
assumption of a field strong enough to raise the elec- 
tron temperature well above that of the heavy particles, 
but that this is usually accomplished even at fairly 
weak field strengths [Ginsburg and Gurevich, I960]. 
Under these conditions the spectral amplitudes do not 
depend explicitly on the strength of the field. For 
elastic collisions, d — 2m/M which is or order 10~ 5 , 
but in an excited gas 8 can become much larger, even 
approaching 10 _1 in some cases (see the reference 
cited above). 

C. Space-charge effects in plasma inhomogeneities: 
Possible quasi-transverse currents at all harmonics. 
The relative amplitudes depend upon both the func- 
tional dependence of no on r and the direction of 
the primary electric field, but are related to the param- 
eter (^m/A w ), where %m is the maximum displacement 
of an electron in the field and A„ is a length measur- 
ing the distance over which the electron density 
changes appreciably. These currents can be expected 
to be most prominent in regions of rapidly changing 
plasma density and in the neighborhood of a point of 
plasma resonance. 



7. References 

\rmstrong, J. \.. Y Bloembergen, .). Ducuing, and P. S. Pershan 
(1962), Interaction between light waves in a nonlinear dielectric 
Phys. Rev. 127, 1918-1940. 
Baird, .1. R., and P. I). Coleman (1961), Frequency conversion in a 

microwave discharge, Proc. [RE 49, No. 12, 1890 1900. 
Bloembergen, N. (1963), Wave propagation in nonlinear electro- 
magnetic media, Proc IEEE 51, No. 1, 124-131. 
Brown, S. C. (1959), Basic Data of Plasma Physics Chap. 1 (John 

Wiley and Sons, New York, N.Y.). 
Chen, K. M. (1962), Interaction of a high-intensity KM field with a 

low-density plasma, IRE Trans. AP-10, No. 1,31-42. 
Epstein, M. (1960), Nonlinear behavior of the electrical conductivity 

of a slightly ionized gas, Phys. Fluids 3, No. 6, 1016-1018. 
Fain, V. M. (1955), The velocity distribution of electrons in the 
presence of a varying electric field and a constant magnetic field, 
J. Exper. Theoret. Phys. USSR 28, 422-430. 
Feinstein, J. (1950), Higher order approximations in ionospheric 

propagation, J. Geophys. Res. 55, 161-170. 
Forsterling, K., and H. Wiister (1951), Uber die Entstehung von 
Oberwellen in der Ionospare, J. Atmospheric Terrest. Phys. 2, 
22-31. 
Ginsburg, V. L. (1958), Nonlinear interaction of radio waves propa- 
gating in a plasma,.!. Exper. Theoret. Phys. USSR 35, 1573-1574. 
Ginsburg, V. L., and A. V. Gurevich (I960). Nonlinear phenomena 
in a plasma located in an alternating electromagnetic Held, Usp. 
Fiz. Nauk 70, 201-246 and 393-428. 
Maiman, T. H. (1957). Solid Stale Millimeter Wave Generation 
Study, Hughes Aircraft Company Physics Laboratory, Final 
Repori on Contract DA 36-O39SC-73063. 
Margenau, H., and L. M. Hartman (1948), Theory of high-frequency 
discharges. II. Harmonic components of the distribution func- 
tion, Phys. Rev. 73, No. 4, 309-315. 
Rosen, 1*. (1961), Generation of the third harmonic by an electro- 
magnetic signal in a plasma, Phys. Fluids 4, No. 3, 341-344. 
Rydbeck, (). F. II. (1963), Dynamic nonlinear electromagnetic wave 
propagation and harmonic radiation in magnetoionic media. Con- 
ference on Wave Interaction and Dynamic Nonlinear Phenomena 
in Plasmas, Pennsylvania Stale University Engineering Proceed- 
ings P-42, 42-70. 
Swan, C. B. (1961), Generation of microwave harmonics in and elec- 
trodeless discharge at low pressure, Proc. [RE 49, No. 12, 1941- 
1942. 
Taylor, L. S. (1961), Conductivity of slightly ionized gases, Phys. 

Fluids 4, No. 12, 1499-1503. 
Uenohara, M., T. Masutani, and K. inada (1957), A new high-power 

frequency multiplier, Proc. [RE 45, 1419-1420. 

Vilenskii, I. \1. (1953), The effect of nonlinearities on radio waves 

propagating in tin* ionosphere, Doklady Akad. Nauk USSR 92, 

525-528. 

Visvanathan, S. (1962), Electron temperature in partially ionized 

gases subject to intense ac fields, J. Appl. Phys. 33, 2481-2483. 

Wet/el, L. (1961a), Wave interaction in plasma inhomogeneities, 

J. Appl. Phys. 32, No. 2, 327-328. 
Wet/el, L. (1961b), Electric-field-induced anisotropics in an inhomo- 

geneous plasma, Phys. Kev. 123, No. 3, 722-726. 
Wetzel, L. (1963), Harmonic excitation in and reradiation from non- 
uniform ionized regions. Symposium on Electromagnetic Theory 
and Antennas, ed. E. C. Jordan, (Pergamon Press, Oxford), 
365-374. 
Whitmer, R. F., and E. B. Barret (1961), Nonlinear interaction of an 
electromagnetic wave with a plasma layer in the presence of a 
static magnetic field. I. Theory of harmonic generation, Phys. 
Rev. 121,661-668. 

(Paper 69D4-494) 



607 



